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Wave turbulence formalism for long internal waves in a 
stratified fluid is developed, based on a natural Hamiltonian 
description. A kinetic equation appropriate for the descrip- 
tion of spectral energy transfer is derived, and its self-similar 
stationary solution corresponding to a direct cascade of energy 
toward the short scales is found. This solution is very close 
to the high wavenumber limit of the Garrett-Munk spectrum 
of long internal waves in the ocean. In fact, a small modi- 
fication of the Garrett-Munk formalism includes a spectrum 
consistent with the one predicted by wave turbulence. 

Introduction. — Remarkably, the internal wave spec- 
trum in the deep ocean has much the same shape wher- 
ever it is observed, except when the observations are 
made close to a strong source of internal waves [1]. This 
observation led Garrett and Munk [2-4] to propose an an- 
alytical form of the internal wave spectrum that approx- 
imates many observations. This spectrum is now called 
the Garrett-Munk (GM) spectra of internal waves. The 
total energy of the internal waves may be represented as 
an integral over spectral energy density 1 



oscillating in the vertical, so that both the hydrostatic 
balance and the WKB approximation apply: 



E = E (k, to) dk dm . 



(1) 



where k and to are the horizontal and vertical compo- 
nents of the wavevector respectively. Under the assump- 
tion of horizontal isotropy, Garrett and Munk proposed 
the following empirical expression for the spectral energy 
density: 



E(k, to) 



2fNEm/m* 



iV) 5/2 {N 2 k 2 + f 2 m 2 ) 



7T 1 



(2) 



Here E is a constant, quantifying the total energy content 
of the internal wave spectrum, / is the Coriolis param- 
eter, N is the buoyancy frequency, k — |k|, and to* is 
a reference vertical wavenumber to be determined from 
observations. 

The dispersion relation underlying this proposed spec- 
trum is that of long internal waves, with a profile rapidly 



: A number of missprints in the original PRL 87 168501 
(2001) manuscript stemed from the fact that GM spectra is a 
one dim spectra, i.e. to obtain total energy one should inte- 
grate it over k and m, not k and m as we mi stakenly thought 
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= f + N 2 (k/ m y 



(3) 



Using this dispersion relation, the spectrum can be trans- 
formed from wave- number space (k,m) into frequency- 
horizontal wavenumber space (k,ui), or frequency- 
vertical wave number space (m,u>). In particular, the 
integral of E(k, to) over fc -or equivalently, the integral of 
E(m, ui) over to- yields the moored spectrum 



E(u) = 2f e(tt (!-(// >u) 2 ) 



1/2 u? 



(4) 



with an 1 /to 2 dependence away from the inertial fre- 
quency that appears prominently in moored observa- 
tions. 

The GM spectrum constitutes an invaluable tool for 
oceanographers, assimilating hundreds of different obser- 
vations into a single, simple formula, that clarifies the dis- 
tribution of the energy contents of internal waves among 
spatial and temporal scales. Yet a number of questions 
regarding the spectrum itself remain open. One is about 
its accuracy: Garrett and Munk basically made up a fam- 
ily of spectra, depending on a few parameters, that was 
simple and consistent with the dispersion relation (3) for 
long internal waves; and then fitted the parameters to 
match observational data. Their success speaks of their 
powerful intuition, yet it leaves the door open to question 
both the accuracy of the parameter fit, and the appro- 
priateness of the proposed family of spectra, necessarily 
incomplete, which involved a high degree of arbitrariness. 
The other question, more fundamental to theorists, is 
to explain or derive the form of the spectrum from first 
principles. Such explanation should surely involve the 
nonlinear interaction among internal modes, as well as 
the nature of the forcing and dissipation acting on the 
system. 

In this work, we elucidate which spectrum the theory of 
wave turbulence (WT) would predict for internal waves in 
scales far away from both the forcing and the dissipation. 
Wave turbulence theory (also called weak turbulence, to 
contrast it to the "strong" turbulence of isotropic fluids) 
applies to Hamiltonian systems characterized by a scale 
separation between a fast, linear dispersive wave struc- 
ture, and its slow, nonlinear modulation. 

In what follows we assume that there is pumping of 
energy into the internal wave field by the wind, by in- 
teraction with surface waves, or by other processes. We 
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assume that these pumping processes can be character- 
ized by wavelengths of the order of hundreds of meters. 
Moreover, we assume that the processes which remove en- 
ergy from the internal wave field, such as wave breaking, 
turbulent mixing, multiple reflections from the surface 
and bottom boundary layers, or interaction with bottom 
topography, can be characterized by lengthscales of less 
than a meter. Then there is a region of lengthscales, 
called inertial interval or transparency region, where dis- 
sipation and pumping processes are not important. It 
is the nonlinear wave interaction which determines the 
form of the spectrum in the region of transparency. En- 
ergy then gets into the system of internal waves at large 
scales, cascades through the inertial scales via multiple 
nonlinear interactions, and is absorbed in the small scales 
by the dissipation processes. This scenario corresponds 
to the Kolmogorov spectrum of WT theory [5] . Accord- 
ing to WT, the energy distribution in the inertial region 
is defined solely by the nature of nonlinear interactions 
and by the linear dispersion relation of the waves in the 
system. 

Notice that, in the transparency region, the dispersion 
relation is self-similar, since \uj\ » f. When this is the 
case, and the nonlinearity is homogeneous, WT theory 
predicts self-similar stationary solutions. For \m\ >> m* 
and \uu\ » /, the Garrett-Munk spectrum (2) becomes 



E(k,m)~ (fc 2 m 3 / 2 ) 



(5) 



is the Lagrangian derivative following a particle. Notice 
that we are considering waves long enough for the hy- 
drostatic balance to be valid, yet not so long to feel the 
effects of the rotation of the earth. This is consistent 
with the scales of the conjectured transparency region 
described above. 

Changing to isopicnal coordinates, where the roles of 
the vertical coordinate z and the density p as indepen- 
dent and dependent variables are reversed, the equations 
become: 

— +u-Vu+ =0, M p = gz, z pt + V-(z p u =0. 

at p 

Here V = (d x ,d y ) is the isopicnal gradient, acting along 
surfaces of constant density, and M is the Montgomery 
potential M = P + pz. We shall consider flows which are 
irrotational along isopicnals; for these, it is convenient 
to introduce a horizontal velocity potential 0, such that 
u = V(/>(x, p, t) , reducing the equations further to the 
pair 



h + j|V*| 2 



+ p/"rn-n<, 

PJ J Pi 

u t + v-(uv<p) =o. (8) 



The spectrum we shall obtain below using the WT for- 
malism is, instead, 



Here we have introduced the variable II = pM pp /g = 
pz p , and, for future convenience, a reference equilibrium 
value n (p) = —g/N 2 , where iV 2 = — gp' {z)/po is the 
square of the buoyancy frequency The variable II, rep- 
resenting the stratification lengthscale, is the canonical 
conjugate of <f) under the Hamiltonian flow given by 



E(k,m)~ (fc 3 / 2 m 3 / 2 ) 



(6) 



The small disparity in the exponents of the Garrett-Munk 
spectrum and the prediction of WT theory may be at- 
tributed either to effects that WT does not capture, such 
as wave interaction with shear and vorticity and wave 
breaking, or to inaccuracies of the GM spectrum itself. 
In the last section of this letter, we introduce a slight 
modification of GM, which yields the spectrum in (6). 

Hamiltonian structure and kinetic equation — The 
equations for long internal waves in an incompressible 
stratified fluid are 

du VP n „ 

-77 + = 0, P z + P g = 0, 

dt p 



dp 



0, V • u + w z 



0. 



where u and w are the horizontal and vertical compo- 
nents of the velocity respectively, P is the pressure, p the 
density, g the gravity constant, V = (d x , d y ) the horizon- 
tal gradient operator, and 

— - — u V w— 

dt dt dz 




n|W| 2 



n-n 



pi 



dp] 



dxdp. (9) 



The first term in this Hamiltonian clearly corresponds 
to the kinetic energy of the flow; that the second term 
is in fact the potential energy follows from the simple 
calculation 



n n 



pi 



dpi 



d P =\ 



(dz — dz ) 



dp 



-(z- z ) 2 dp - 



-p (z -z )dz + d[ ^p (z - z ) 2 ) , (10) 



so that the second term in (9) is simply g j p(z — Zo)dz . 
The equations of motion (7, 8) can be written in terms 
of the Hamiltonian (9) in the canonical form 



m. = 



SH 



SH 



(ii) 



For simplicity, we shall take the buoyancy frequency 
TV of the equilibrium profile to be a constant, and we 
shall replace the density p in the denominator of the 
Hamiltonian's potential energy by a constant p . This is 
the WKB -also Boussinesq- approximation, which makes 
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sense for waves varying rapidly in the vertical direction, 
particularly since the water density typically changes 
only by a few percent over the full depth of the ocean. 

Let us decompose II into the sum of its equilibrium 
value and deviation from equilibrium II = n +II' . Since, 
in the WKB limit, the linear part of the resulting Hamil- 
tonian has constant coefficients, it is natural to perform 
a Fourier transformation in both vertical and horizon- 
tal directions. We assume that both Fourier spectra are 
continuous, which is a reasonable approximation if the 
wave's vertical wavelengths are much smaller than the 
depth of the ocean. Then 



H 



27V 2 ' 



9 in' i 2 

2p 2 m 2 p 



dp 



(2tt) 3 / 2 



k 2 • k 3 hp^pj^ 



<KPl + P2 - P3) rfpi23 , 



(12) 



where p is a three-dimensional wave vector p = {k, m} 
and dpi23 = rfpirfp2^P3- 

We now introduce a canonical transformation which 
transfers the equations of motion for the Fourier compo- 
nents (f> p and lip into a single equation for the canonical 
variable a p . We choose this transformation so that the 
quadratic part of Hamiltonian becomes diagonal in a p . 
The canonical transformation reads: 



N i * \ tt/ Pomk t 

1 1 . r(a P - o_ p ), n p = \/^-(a p + a_ p ), 



2p mk 



Then the pair of the canonical equations of motion (11) 
become the single equation 



.d_ _&H_ 
l dt Up ~ da* p ' 



(13) 



with Hamiltonian 



H = J uj p \a p \ 2 dp + 
J Vp iP2 p 8 (a* pi a p2 a p3 + a pi a p2 a* p3 ) 5 pi+p2 _ p3 dp 123 + 

j ^PlP2P3 ( a pi a P2 a p3 + a Pl fl P2 a P 3 ) <5pi+P2+P3 ^Pl23- 



Here u> p is the linear dispersion relation for the Hamilto- 
nian (12), 

9 k 



and Vp lP2 p 3 
ment, given by 



Np \m\ 

is the internal wave interaction matrix ele- 




The form of this Hamiltonian is typical for systems 
with three-wave interactions and cylindrical symmetry. 
Following wave turbulence theory, one proposes a per- 
turbation expansion in the amplitude of the nonlinear- 
ity To zeroth order in the perturbation, one recovers 
the linear waves. At higher orders, the nonlinear interac- 
tions lead to a slow modulation of the wave amplitudes, 
representing spectral transfer of the conserved quanti- 
ties of the Hamiltonian. This transfer manifests itself in 
the perturbation expansion through resonances or secu- 
lar terms, occurring on the so called resonant manifold. 
Energy transfer is described by an approximate kinetic 
equation for the "number of waves" or wave- action n p , 
defined by 

n p S(p - p') = (a*ap/) . 

This kinetic equation is the classical analog of the Boltz- 
mann collision integral; it has been used for describibg 
surface water waves since pioneering works by Hassel- 
mann [6] and Zakharov [7,8]. The derivation of the ki- 
netic equation using the wave turbulence formalism can 
be found, for instance, in [5]. For the three- wave Hamil- 
tonian (14), the kinetic equation reads: 



dn 
~~dt 



PP1P2 1 



/p!2 S p - 



8u B 



-2ir J |Vp iP p 2 | 2 /ip2<5pi-p- 



P2 u w P1 — cj p — cjp 2 



,dpi2 



dPl2 , (15) 



where f pl2 = n pi n P2 - n p (n pi + n P2 ) . 

Assuming horizontal isotropy, one can average (15) 
over all horizontal angles, obtaining 

dn If 

= — / [R\2 ~~ R\i ~~ -^lfe) dkidk2dmidm,2 , 
dt k J 

R i2 = A ki2 5 ( u p - u pi - ^pi)I\2 I V12 1 2 S m -m 1 - m2 kk 1 k 2 , 
A ^ = (<5(k ki k a )> = f S(k - ki - k a ) de 1 de 2 , 

Afci2 = ^2(^)2 + (kk 2 ) 2 + (fcife) 2 ) -k A -k\ — k\. (16) 

In wave turbulence theory, three-wave kinetic equa- 
(l4) ons admit two classes of exact stationary solutions: 
thermodynamic equilibrium and Kolmogorov flux solu- 
tions, with the latter corresponding to a direct cascade of 
energy -or other conserved quantities- toward the higher 
modes. The fact that the thermodynamic equilibrium -or 
equipartition of energy- n p — 1/uj p is a stationary solu- 
tion of (16) can be seen by inspection, whereas in order 
to find Kolmogorov spectra one needs to be more elab- 
orate. Let us assume that n p is given by the power-law 
anisotropic distribution 



rik,m = k x \m\ y . 



(17) 



We will find exponents x and y by requiring that (17) is 
a stationary solution to (16). We shall use a version of 
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Zakharov's transformation [7,8] introduced for cylindri- 
cally symmetrical systems by Kuznetsov in [9]. Let us 
subject the integration variables in the second term R\ 2 
in (16) to the following transformation: 

fci = k 2 /k[, mi — m 2 /m'i, k<i — kk' 2 /k'i, m,2 — mm' 2 /'m l 1 . 

Then R\ 2 becomes i?f 2 multiplied by a factor 



-6-2z / \ 2+2y 

m \ 
mi) 



Furthermore, let us subject the third term R 2 k in (16) to 
the following conformal transformation: 

ki = kk[/k' 2 , mi — mm[/m 2 , k 2 — k 2 /k' 2 , m 2 = m 2 /m' 2 . 

Then R\ k becomes R\ 2 multiplied by a factor 



-6-2a; 



(-) 



2+2y 



Therefore (16) can be written as 



dt k 



mi 



, \ -6-2x / s 2+2^ 

k 2 \ I m 

m 2 



dkidk2<imidm2 . (18) 



We see that the choice — 6 — 2x = 2 + 2y = 1, which 
gives x = — 7/2, y = — 1/2 , makes the right-hand side 
of (16) vanish due to the delta function in the frequen- 
cies, corresponding to energy conservation. The resulting 
wave action and spectral energy distributions are given 
by 



Ikl-^lmr 1 / 2 , 



(19) 



|k| 



- 3 / 2 lm|- 3 / 2 , 



This solution, ((6) of the introduction) corresponds to 
the flux of energy from the large to the small scales. 

This is the main result of this article - the short wave 
part of the Garrett-Munk spectra is close to the station- 
ary solution to a kinetic equation hereby derived for inter- 
nal waves, based on a Hamiltonian structure appearing 
naturally in isopicnal coordinates. 

A modified Garrett-Munk spectrum — In this section, 
we show how a minor modification of the Munk-Garrett 
formalism allows us to match the wave turbulence pre- 
diction (6). As in [3], we introduce two functions 



A(X) 



t-1 



B(co) 



2/ 



1 



- (i)" 



1/2 



and a reference, frequency dependent wave vector 
(k*,m*). In [3], m* was a constant, and k* was given 



by k* = (w 2 — f 2 ) 1 / 2 ™*. However, the same formalism 
carries through if one allows both k* and m* to depend 
on uj, provided that the condition k*/m* = \J w 2 — f 2 is 
met. In keeping with the spirit of self-similarity, we shall 
propose that 



1-5 

2 j-2\-2~ 



fc*= 7 (c 2 -/ 2 ) 



m 



7 f) 1 , (20) 



where 5, like t and 7, is a constant to be determined 
from observations (or, in our case, from wave turbulence 
theory). 

Then, following [3], we propose an energy spectrum of 
the form 



E(k,u>) = EA 



k \ B(lu) 
k* 



It follows from the dispersion relation that we have 



B(u) 



E{k,m) = 



y ' ' \m*J m 

2fNE {m/m*)A{m/m* 



N 2 k 2 + f 2 m 2 



(21) 
(22) 



If we pick t = 2 and S = —1/2, the asymptotic behavior 
of (22) agrees with the prediction (6) of wave turbulence 
theory. In [3], 5 was zero by default, and t = 2.5. Notice 
that, independently of the choices of t and d, the moored 
spectrum is always given by (4). 

Conclusion. — We have found a natural Hamiltonian 
formulation for long internal waves, and used it within 
the wave turbulence formalism to determine the station- 
ary energy spectrum corresponding to a direct cascade 
of energy from the long to the short waves. This spec- 
trum is close to the one that Garrett and Munk fitted to 
available observational data. The small difference could 
be due either to physical effects that the wave turbulence 
formalism fails to capture, or to a real necessary correc- 
tion to the GM spectrum. We show how a slight modi- 
fication of the GM spectrum yields results in agreement 
with WT theory. 
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